Lowest Order Constrained Variational 



Calculation of the Polarized Nuclear Matter 
with the Modern AVis Potential 

G.H. Bordbar *tand M. Bigdeli 
Department of Physics, Shiraz University, Shiraz 71454, Iran^ 

and 

Research Institute for Astronomy and Astrophysics of Maragha, 
P.O. Box 55134-441, Maragha, Iran 

Abstract 

The lowest order constrained variational method is applied to calculate the po- 
larized symmetrical nuclear matter properties with the modern AVis potential per- 
forming microscopic calculations. Results based on the consideration of magnetic 
properties show no sign of phase transition to a ferromagnetic phase. 
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1 Introduction 

The properties of dense matter is a subject that theoretical physicists have desired to 
study. The magnetic property of nucleon matter is of special interest in nuclear and 
astrophysics which can be related directly with magnetic source of pulsars, rapidly rotat- 
ing neutron stars with strong surface magnetic fields in the range of 10^^ — 10^^ Gauss 
[1, 2, 3]. The most interesting and stimulating mechanisms that have been suggested 
is the possible existence of a phase transition to a ferromagnetic state at densities cor- 
responding to the theoretically stable neutron stars and, therefore, of a ferromagnetic 
core in the liquid interior of such compact objects. Such a possibility has been studied 
by several authors using different theoretical approaches [4-24], but the results are still 
contradictory. Vidana et al. [21], Vidana and Bombaci [22] have considered properties 
of spin polarized neutron matter and polarized isospin asymmetric nuclear matter using 
the Brueckner-Hartree-Fock (BHF) approximation by employing three realistic nucleon- 
nucleon interactions, Nijmegen II, Reid93 and NSC97e respectively. Zuo et al. [25] have 
also obtained properties of spin polarized neutron and symmetric nuclear matter using 
same method with AVis potential. The results of those calculations show no indication 
of ferromagnetic transition at any density for neutron and asymmetrical nuclear matter. 
Fantoni et al. [20] have calculated spin susceptibility of neutron matter using the Auxiliary 
Field Diffusion Monte Carlo (AFDMC) method employing the AU6 + UIX three-body 
potential, and have found that the magnetic susceptibility of neutron matter shows a 
strong reduction of about a factor 3 with respect to its Fermi gas value. Baldo et al. [26] , 
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Akmal et al. [27] and Engvik et al. [28] have considered properties of neutron matter 
with AVis potential using BHF approximation both for continuous choice (BHFC) and 
standard choice (BHFG), variational chain summation (VCS) method and lowest order 
Brueckner (LOB) respectively. On the other hand some calculations, like for instance the 
ones based on Skyrmehke interactions predict the transition to occur at densities in the 
range (1 — 4)po (po = 0.16/m~^) [29] . This transition could have important consequences 
for the evolution of a protoneutron star, in particular for the spin correlations in the 
medium which do strongly affect the neutrino cross section and the neutrino mean free 
path inside the star [30]. 

Recently, we have used the lowest order constrained variational (LOCV) method [31] to 
calculate the equation of state of symmetrical and asymmetrical nuclear matter and some 
of their properties such as symmetry energy, pressure, etc. [32-35]. We have also obtained 
the properties of spin polarized liquid ^He [36] using this method. The LOCV method is 
a useful tool for the determination of the properties of neutron, nuclear and asymmetric 
nuclear matter at zero and finite temperature. It is a fully self-consistent formalism which 
does not bring any free parameters into calculation. It employs a normalization constraint 
to keep the higher order term as small as possible [31]. The functional minimization pro- 
cedure represents an enormous computational simplification over unconstrained methods 
that attempt to go beyond lowest order. 

In our pervious work, we have developed the LOCV method to compute the properties 
of polarized neutron matter such as total energy, magnetic susceptibility, pressure, etc. 
[37] , and have seen that the spontaneous phase transition to a ferromagnetic state in the 
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neutron matter does not occur. In the present work, we intend to calculate the polarized 
symmetrical nuclear matter properties using the LOCV method with the modern AVis 
potential [38] employing microscopic calculations. 



2 LOCV FORMALISM 

We consider a cluster expansion of the energy functional up to the two-body term, 

The smallness of the three-body cluster energy has been discussed in Ref . [32] , where it 
is shown that our cluster expansion converges reasonably and it is good approximation 
to stop after the two-body energy term. This property can also be predicted by looking 
at the correlation between the particles which will be discussed in the next section. 

The one-body term Ei can be written as Fermi momentum functional {kp = (37r^p^*'')^), 

Labels 1 and 2 are used instead of spin up and spin down nucleons, respectively, and 
p — p^^p^\s, the total nuclear matter density. The two-body energy E2 is 

^2 = -^-E(^J>(12)|U-J^), (3) 



2^ 



where 



""^^^^ ^ "i^^t^^^^^' t^'^' ^^^^^^^ + /(12)l-(12)/(12) (4) 



/(12) and V{12) are the two-body correlation and potential, respectively. For the two- 
body correlation function, /(12), we consider the following form [32, 33]: 



/(12) = ^/(^)(12)0('=)(12), 



(5) 



k=l 



where, the operators 0(*^)(12) are given by 



(6) 



and 5'i2 is the tensor operator. A complete discussion of correlation function and especially 
its form are given in Ref. [31]. 

After doing some algebra, we find the following equation for the two-body energy of 
the polarized symmetrical nuclear matter. 



Eo — 



-0'zl-0'z2 I SSz 



(2J+1)(2T + 1) ^^^s+^ 

^'P\'^^) JLTSS. 2(25+1) ^ ^ ^ J 

xjdr[ [f^^"a^f{kfT) + ^ ({K - + (K - 3K.)(4r - 3) 

(1)2. 
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+ E /i'^ +^({K + K + (-6A; + 14)T4 + -(A;-l)l^, 
+[K, + l^,, + (-6A; + U)Vtz - (A) - l)Vi,,](4r - 3) 
+ [Vt + + (-6A; + 14)V^iT][4T]} a«'(V) 

+ [Vi2 + 1^2. + (V^2r + 1^2.r)(4r - ^)]cf{kfr) 



+ [(^.2 + W(4r-3)]ci(f) (V) /i' 



2777, 

+ ^[[(^- - 2(V^/2.r + Vl2r) " 3V^,2r)(4T - 3)] 



{J, L, 5", Sz} and the coefficient a^^^^, etc. are defined as follows, 



^ 2J+1' ' 2J+1' 2J+1 



^i = L(L + l); ^. = ^j^; ^3 = (18) 



J(J2 + 2J+1) j(j2 + j + 2) 

2J + 1 ' 2J + 1 (''^ 



J^ + 2J^ + 2J+1 _ J(J^ + J + 4) 
2JTi ' 2J+1 



2J + 1 ' 2JTi ^^^^ 



and 



Ij,sM = J dqPsM)JjM (22) 
In the above equation, Jj{x) is the Bessel's function and Ps^{q) is defined as follows: 

-^((^^"r-(^F^^)')V^ + g'] (23) 
for I - A;^^^| < g < i + k''p?'\, 

PsM-l^rnin{kr,k^F') (24) 

for q < i|A;J-i ^^^^^ 

i'5.(g)=0 (25) 

for g > i + fc^^'^l, where a^i and (T22 are equal to |, — | for spin up and spin down 
nucleons, respectively. 
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Now, we minimize the two-body energy, Eq.(7), with respect to the variations in the 
correlation functions fa^''\ but subject to the normahzation constraint [33], 



jUij\hl-f{l2)\ij)a^0, 



(26) 



where in the case of spin polarized nuclear matter, the function hs^ (r) is defined as 

21 -1/2 



h.(r) 



±1 



(27) 



1 ^. = 

Prom the minimization of the two-body cluster energy, we get a set of coupled and 
uncoupled differential equations which are the same as presented in Ref. 



3 RESULTS 

In Fig. 1, we have shown the correlation function versus the relative distance (r). Fig. 
1 shows that the correlation between particles is short range and heals to 1 very quickly. 
This means that the two-body term mainly contributes to the interaction of particles and 
therefore higher order terms can be neglected. 

The energy per particle of the polarized symmetrical nuclear matter versus density 
for different values of the spin polarization have been shown in Fig. 2. This figure shows 
that the low polarization gives more binding energy than the high polarization. It is also 
seen that there is no crossing of the energy curves of different polarizations, vice versa 
by increasing density, the difference between the energy of nuclear matter at different 
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polarization becomes more sizable. This shows that the spontaneous phase transition to 
a ferromagnetic state in the symmetrical nuclear matter does not occur. 

In Fig. 3, we have plotted the quadratic spin polarization dependence of energy 
per particle at different densities. As can be seen from this figure, there are two points 
worth stressing. First the energy per particle of the polarized symmetrical nuclear matter 
increases as the polarization increases and the minimum value of energy occurs at 5 — 
for all densities. This indicates that the ground state of symmetrical nuclear matter is 
paramagnetic. Second the variation of the energy of symmetrical nuclear matter versus 
is nearly linear, 



In Fig. 3, the results of ZLS calculations [25] are also given for comparison. There is an 
agreement between our results and those of ZLS, specially at low densities. 

The magnetic susceptibility, x, which characterizes the response of a system to the 
magnetic field and gives a measure of the energy required to produce a net spin alignment 
in the direction of the magnetic field, is defined by 



where M is the magnetization of the system per unit volume and H is the magnetic 
field. We have calculated the magnetic susceptibility of the polarized symmetrical nuclear 
matter in the ratio x/Xf form. By using the Eq. 29 and some simplification, the ratio of 



E{p, 5) 



E{p,0) + anuci{p)5'^- 



(28) 




(29) 
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X to the magnetic susceptibility for a degenerate free Fermi gas xf can be written as 

where Ep = 'h?k\/2m is the Fermi energy and kp = (3/27r^p)^/^ is Fermi momentum. 
Our results for magnetic susceptibility are displayed as a function of density in Fig. 4. As 
can be seen from Fig. 4, this ratio changes continuously for all densities and decreases as 
the density increases. Therefore, the ferromagnetic phase transition is not predicted by 
our calculation. For comparison, we have also shown the results of ZLS [25] in this figure 
which shows good agrement whit our results. 

By differentiating symmetrical nuclear matter energy curve at each polarization [5) 
whit respect to the density we can evaluate the corresponding pressure, 

P(p,.)=p^^, (31) 

In Fig. 5, we have shown the pressure of polarized symmetrical nuclear matter as a 
function of density p for various polarizations. We see that equation of state of polarized 
symmetrical nuclear matter, P(p, 5), becomes stiffer by increasing the polarization in the 
density range which was considered. 

In Fig. 6, we have also presented the Landau parameter, Go, which describes the spin 
density fluctuation in the effective interaction, versus density. It is seen that the value 
of Go is always positive and monotonically increasing up to highest density and does 
not show any magnetic instability for the neutron matter. A magnetic instability would 
require Go < —1. 
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4 Summary and Conclusions 

We have computed the magnetic properties of polarized symmetrical nuclear matter, that 
is related directly with magnetic source of pulsars, and other properties using the lowest 
order constrained variational (LOCV) method with with AVig potential. We have stud- 
ied the total energy per particle of nuclear matter as a function of density and the spin 
polarizations 6. We have found that in the range of densities explored, difference between 
the energy of polarized nuclear matter at different polarization becomes more appreciable. 
We have also seen that total energy per particle is parabolic on the spin polarization 5 in a 
very good approximation up to full polarization for all densities. Magnetic susceptibility, 
which characterizes the response of the system to the magnetic field was calculated for the 
system under consideration and was found that it changes continuously for all densities. 
There is an overall agreement between our result and those of Zuo et al. [25] . In conclu- 
sion, we see that equation of state of polarized symmetrical nuclear matter becomes stiffer 
by increasing the polarization in the density range which was considered. The Landau 
parameter, Gq has been considered and it was seen that the value of Go is always positive 
and monotonically increasing up to high densities. Finally, our results have shown no 
phase transition to ferromagnetic state. 
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Figure 1: The two-body correlation functions of the full polarized nuclear matter as a 
function of relative distance at p = 0.67/m~^. 
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Figure 2: The energy per particle of the polarized symmetrical nuclear matter versus 
density(p) for different values of the spin polarization (6). 
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Figure 3: Our results (full curves) for the energy difference of polarized and unpolarized 
cases versus quadratic spin polarization (5^) for different values of the density(p) of the 
neutron matter. The results of ZLS [25] (dashed curves) are also presented for comparison. 
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Figure 4: Our result (full curve) for the magnetic susceptibility of the polarized symmetri- 
cal nuclear matter as the function of density(p). The results of ZLS [25] (dashed curves) 
are also given for comparison. 
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Figure 5: The equation of state of polarized symmetrical nuclear matter for different 
values of the spin polarization (5). 
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Figure 6: Our result for the Landau parameter, Go, as function of density(p). 
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